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ABSTRACT 

Spectra of the cellular photospheric flows are determined from full-disk Doppler velocity 
observations acquired by the Helioseismic and Magnetic Imager (HMI) instrument on the Solar 
Dynamics Observatory (SDO) spacecraft. Three different analysis methods are used to separately 
determine spectral coefficients representing the poloidal flows, the toroidal flows, and the radial 
flows. The amplitudes of these spectral coefficients are constrained by simulated data analyzed 
with the same procedures as the HMI data. We find that the total velocity spectrum rises 
smoothly to a peak at a wavenumber of about 120 (wavelength of about 35 Mm), which is typical 
of supergranules. The spectrum levels off out to wavenumbers of about 400, and then rises again 
to a peak at a wavenumber of about 3500 (wavelength of about 1200 km), which is typical of 
granules. The velocity spectrum is dominated by the poloidal flow component (horizontal flows 
with divergence but no curl) at wavenumbers above 30. The toroidal flow component (horizontal 
flows with curl but no divergence) dominates at wavenumbers less than 30. The radial flow 
velocity is only about 3% of the total flow velocity at the lowest wavenumbers, but increases in 
strength to become about 50% at wavenumbers near 4000. The spectrum compares well with 
the spectrum of giant cell flows at the lowest wavenumbers and with the spectrum of granulation 
from a 3D radiative-hydrodynamic simulation at the highest wavenumbers. 

Subject headings: convection, Sun: granulation, Sun: interior, Sun: photosphere 

1. INTRODUCTION sphere. While buoyancy and heat transport drive 

these motions, these flows also transport magnetic 
The Sun is the only star for which we can di- field anc j an g U i ar momentum and thus play a crit- 

rectly observe the convective motions that carry ical ro i e j n ^g ro tational and magnetic history of 

energy fiom its ladiative interior to its photo- the star. The convective motions, and their in- 
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teractions with magnetic fields, are the ultimate 
drivers behind many scientific questions concern¬ 
ing the Sun what is the source of the 11- 
year sunspot cycle, what heats the Sun’s corona, 
what accelerates the solar wind, what triggers so¬ 
lar flares, prominence eruptions, and coronal mass 
ejections? 


Solar granules, bright cells with darker borders, 
were observed with solar telescopes by the 1860s. 
These cellular flows have diameters of about IMm, 
flow speeds of about 3000 m s -1 , and lifetimes of 
about 10 minutes ( Brav et al]|l984j ). 

The much l arger supergranules were first no¬ 
ticed bv lHartj f 19541) but named and better char¬ 
acterized by Leightoiretjd] ( 1963) • Supergranules 
have diameters of about 35 Mm, flow speeds of 
abou t 40 0 m s" 1 , and lifetimes of about a day 
( Rieutord fc Rinconl 2010 ). 

Even larger cells, giant cells, were pro posed to 
exist based on theoretical arguments (|Simon & Weis s 


19681) . Hints of the possible existence of these cells 


have been uncovered over the ensuing decades but 
detailed observatio ns and characterizatio ns have 
only come recently ( Hathaway et al. 2013t h These 
cells have widths of 100 Mm or more, flow speeds 
of about 10 m s -1 , and lifetimes of several months. 


Yet an other set of cel ls , mes ogranules, were 
found by iNovember et al.1 ( 198ll ) to be interme¬ 
diate in size, flow speeds, and lifetimes between 
granules and supergranules. 

The advent of space based observations, with 
both continuous coverage and high spatial and 
temporal resolution, has enabled many new 
studies of solar convection. Data from the 
Michelson Doppler Investigation (MDI) on the 
ESA/NASA Solar and Heliospheric Observatory 
(SOHO) mission were analyzed by several groups. 
lLawrence et al. (1999;) used a wavelet analysis on 
Doppler, intensity, and magnetic images. They 
noted that granules dominated the radial flows 
while supergranules dominated the horizontal 
flows, but mesogranules were notably weak. (They 
also concluded that their wavelet analysis only 
captured about 30% of the signal.) The MDI 
Doppl er data were analyzed and simulated by 
Hathaway et al. ( 2000h to obtain photospheric ki¬ 


netic energy spectra covering spherical harmonic 
wavenumbers, £, from i ~ 1 to l ~ 3000. They 
concluded that only two spectral features were 


evident, a peak representing supergranules at 
i ~ 120 and a broad pe ak representing granu les at 
t > 1000. Others (e.g. iRieutord et al. l2000h have 
also concluded that mesogranulation is not a true 
or discrete scale of solar convection. 


Power spectra of th e flows at very high wav enum- 
bers were obtained by Rieutord et al . ( 20ld l from 
high resolution images obtained with the Solar 
Optical Telescope (SOT) on the Hinode space¬ 
craft. Their horizontal flow spectra show a promi¬ 
nent peak in power at wavenumbers typical for su¬ 
pergranules, along with a continuum of power out 
to higher wavenumbers. (Their correlation track¬ 
ing method for measuring the horizontal flows lim¬ 
ited their measurements to cells larger than ~ 2.5 
Mm or l ~ 1800.) Their vertical flow spectra 
(obtained from Doppler velocity measurements) 
show a prominent peak in power at wavenumbers 
typical of granules {t ~ 2500) with a rapid fall off 
at higher wavenumbers. 


Here we investigate the Sun’s photospheric con¬ 
vection spectrum by analyzing Doppler data from 
the HMI instrument. These data have much 
higher spatial resolution, and better instrumental 
characterization, than the MDI data. The higher 
resolution data gives better confidence in char¬ 
acterizations of the high-wavenumber flow struc¬ 
tures. The better characterization of the HMI 
instrument, its imaging artifacts, and its point- 
spread-function, gives better confidence in the 
characterizations at both ends of the spectrum. 


We employ a new Doppler data analysis tool 
to determine the toroidal component of the hori¬ 
zontal flows (vortical flows with curl but no di¬ 
vergence). While the short (compared to the 
27-day solar rotation period) lifetimes of su¬ 
pergranules suggest that solar rotation should 
not be a big influence on th e flows, numerical 
simulations ( Hathaway! 119821) and observations 
( Duvall fc Gizon 2000h indicate that the flows in 
supergranules have a preferred sense of vortic- 
ity in each hemisphere. This tendency should be 
stronger for the larger, slower cellular flows. 


In the following section we describe the data 
and the data preparation. These data directly 
measure only the line-of-sight velocity component 
and only over the visible hemisphere. We analyze 
the data with three different analysis methods in 
an effort to separate the Doppler signal into three 
components associated with the full vector veloc- 
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ities in the convective flows. In order to fully, 
and accurately, quantify the vector velocities we 
construct simulated data in which we know the 
actual vector velocities over the entire solar sur¬ 
face. The simulated data are then run through the 
same analysis procedures, results are compared 
with those from the HMI data, and the simulated 
data are then revised until a match is obtained. 


This data simulation is described in Section 3. 
The Doppler velocity spherical harmonic spectrum 
analysis is described in Section 4. An analysis that 
separates the radial flow component from the hor¬ 
izontal flow component is described in Section 5. 
An analysis that separates poloidal flows (diverg¬ 
ing horizontal flows) from toroidal flows (horizon¬ 
tal flow vortices) is described in Section 6. The full 
spectrum of motions for all flow components is pre¬ 
sented in the concluding section along with com¬ 
parisons to the spectrum from a numerical simula¬ 
tion of the smallest length scales, and to the spec¬ 
trum of giant cell flow velocities from the corre la¬ 


tion tr acking method described in lHathawav et al 

jaoia . 


2. DATA PREPARATION 


The data used in this study consist of line-of- 
sight Doppl er velocities measure d by the HMI 
instrument ( Scherrer et al.1 12012 1 on the SDO 
spacecraft. These 4096-by-4096 pixel images 
are obtained every 45 seconds and then aver¬ 
aged over 720 seconds t o remove th e 5-m inute 
p-mode oscillati o n sign al. ISchou et al.1 (1201211 and 


ICouvidat et al. ( 2012 ) describe the data acqui¬ 
sition and calibration procedures. The Doppler 
velocity is calculated from the intensity measured 
at six different spectral tuning positions around 
the wavelength of the selected Fe I absorption line 
at 6173 A. An example of a data image is shown 
in Figure [U 


Our first step is to remove the velocity signal 
due to the motion of the spacecraft and the veloc¬ 
ity artifacts associated with the instrument itself. 
Using the spacecraft velocity relative to the Sun 
still leaves behind diurnal signals clearly associ¬ 
ated with the spacecraft’s geosynchronous orbit. 
We found that these variations alter the artifact 
signal as a function of spacecraft velocity toward 
or away from the Sun. The artifacts also change 
when the instrument focus is changed. To rnini- 



Fig. 1.— Doppler velocity image from the 
SDO/HMI instrument for November 1, 2012 at 
5:00UT. The line-of-sight velocity signal ranges 
from -3000 m s -1 (dark blue) to +3000 m s -1 
(dark red). The signal is dominated by the ± 
2000 m s -1 signal from the Sun’s rotation but the 
mottled pattern produced by the solar convection 
is also very apparent (as are “blemishes” due to 
flows around sunspots). 


mize the impact of these variations on our results, 
we chose data taken when the spacecraft motion 
toward or away from the Sun was less than 300 m 
s _1 (one tenth of the ± 3000 m s -1 daily range). 

We extract the imaging artifacts by first remov¬ 
ing the large-scale solar velocity signals that are 
fixed in the image (e.g. the Doppler signal due to 
the axisymmetric flows, the gravitational red shift, 
and the convective blue shift). We then average 
over many hundreds of images so that the signal 
due to the more or less random convective flows 
cancel as they evolve and rotate across the visible 
disk while the artifacts remain fixed. The arti¬ 
fact image obtained from data acquired between 
the dates of October 1, 2012 and May 31, 2013 
is shown in Figure [2] These artifacts would intro¬ 
duce spurious features at low spatial wavenumbers 
in our spectral analyses. 

For this study we have analyzed data acquired 
in November 2012. After removing the Doppler 
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Fig. 2. Doppler velocity image showing imaging 
artifacts from the SDO/HMI instrument for Oc¬ 
tober 2012 through May 2013. The line-of-sight 
velocity signal ranges from -150 m s -1 (dark blue) 
to +150 m s -1 (dark red). 


signals due to the spacecraft motion and the imag¬ 
ining artifacts, each Doppler im age is an alyz ed us- 
ing the procedures de scribed in lHathawavl (1198711 
and lHathawavl ( 1992h to measure and remove the 
Doppler signals due to the axisymmetric flows (dif¬ 
ferential rotation and the meridional circulation) 
and the signal due to the convective blue shift. 

The convective blue shift is an apparent flow 
toward the observer caused by the correlation be¬ 
tween radial upflows and emergent radiative inten¬ 
sity in the unresolved convective flows (granules). 
(The Doppler velocity is an intensity weighted av¬ 
erage over the flow structures within a pixel.) It 
varies systematically from disk center to limb as 
the line-of-sight samples different velocity compo¬ 
nents and different levels of the solar atmosphere. 
The convective blue shift changes substantially in 
regions where magnetic fields thread through the 
surface. Therefore, the convective blue shift sig¬ 
nal requires additional corrections. Figure [3] shows 
the convective blue shift signal for both field free 
regions and for regions with varying magnetic flux 
density. (Note that while a 636.4 m s -1 gravi¬ 
tational red shift has been subtracted from the 


signal, the original zero point was arbitrarily set 
to make the data median match the spacecraft’s 
radial velocity.) The convective blue shift is sub¬ 
stantially altered in magnetic areas. Uncorrected, 
this introduces Doppler signals that are not di¬ 
rectly related to actual flows on the Sun’s surface. 
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Fig. 3.- The convective blue shift signal from the 
SDO/HMI instrument for regions free of magnetic 
field (black line) and regions with radial magnetic 
flux density from 100 to 200 G (blue line), 200 to 
400 G (red line), and 400 to 1200 G (green line). 

We fit the convective blue shift signal, CBS(x), 
to shifted Legendre polynomials (polynomials or¬ 
thonormal with uniform weights on the interval 
0 < x < 1) with the dependent variable 

x = 1 — cos p (1) 

where p is the heliocentric angle of a point from 
disk center. The first five shifted Legendre poly¬ 
nomials are: 


P*(x) = 1 (2) 

P*{ x) = y/Z(2x-l) (3) 

P*( x ) = y/l(<ox 2 - 6x + 1) (4) 

P*(x) = v / 7(20x 3 - 30a; 2 + 12a; - 1) (5) 

P 4 * (x) = 3(70x 4 - 140x 3 + 90x 2 - 20x + 1) (6) 
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We find that the convective blue shift signal is 
well fit with these first five polynomials such that 


example of these corrected data is shown in Figure 
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CBS(x) = Y, C nP:(x) ( 7 ) 

n —0 

The fit coefficients for the different magnetic 
flux densities are given in Table [I] 


Flux Density 

Go 

c\ 

c 2 

Gs 

g 4 

< 100 

-584 

105 

92 

19 

0 

[100, 200] 

-508 

35 

113 

37 

3 

[200,400] 

-490 

-4 

123 

59 

15 

[400,1200] 

-497 

-35 

102 

71 

33 


Table 1: Convective blue shift fit coefficients (in 
units of m s -1 ) for increasing magnetic flux den¬ 
sity (in units of Gauss). 


We correct the data for the convective blue shift 
differences in magnetic regions and find that a sin¬ 
gle convective blue shift signal, shown in Figured 
represents all the data. We then subtract this sig¬ 
nal from the data. 
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Fig. 5.— Prepared Doppler velocity image from 
the SDO/HMI instrument for November 1, 2012 
at 5:00UT. The line-of-sight velocity signal ranges 
from -1000 m s -1 (dark blue) to +1000 m s -1 
(dark red). 


3. THE DATA SIMULATION 


We produce simulated data (for which we know 
the full vector velocities over the entire surface 
of the Sun) that faithfully reproduce the HMI 
Doppler data. 


We represent the vector velocity field on the 
surface of a sph ere by a spectrum of poloidal and 
toroidal modes ( Chandrasekharlll96ir i with 


Fig. 4.— The convective blue shift signal after cor¬ 
recting for magnetic effects. The four lines repre¬ 
sent regions free of magnetic field (black line) and 
regions with radial magnetic flux density from 100 
to 200 G (blue line), 200 to 400 G (red line), and 
400 to 1200 G (green line). 

After removing these signals (the convective 
blue shift, the artifacts, the Doppler signal due 
to differential rotation and meridional flow) the 
Doppler data are dominated by the signal from 
the line-of-sight velocities in the cellular flows. A 
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where Y™(0,(f>) is a spherical harmonic function 
of degree t and azimuthal order m, 9 is the colati¬ 
tude measured southward from the north pole, and 
</> is the azimuth measured prograde from the cen¬ 
tral meridian. The complex coefficients R™, S™, 
and T™ are the spectral coefficients for the radial, 
poloidal, and toroidal components, respectively. 

We simulate the observed line-of-siglit veloc¬ 
ity by specifying these complex spectral coeffi¬ 
cients for all spectral components up to a max¬ 
imum spherical harmonic degree, £ ma ,x = 4096. 
This is an iterative, trial and error, process in 
which spectral amplitudes are prescribed, simu¬ 
lated data are constructed and analyzed, the re¬ 
sults are compared to those from the HMI data, 
and new spectral coefficients are constructed based 
on those comparisons. We start with an estimate 
of the spectral amplitudes as functions of l. The 
complex coefficients for each azimuthal order, to, 
use these spectral amplitudes but with random 
phases for the real and imaginary parts. We use 
the same phase for the radial and poloidal spectral 
coefficients since these two components are linked 
through the mass continuity equation (V • v = 0) 
with 


e(e + i)sr = ~[r 2 R?(r)] (11) 

We do not link the phase of the toroidal compo¬ 
nent to the other two. 

We calculate the three vector velocity compo¬ 
nents from this spectrum at an array of points 
with 4096 points equi-spaced in 9 and 8192 points 
equi-spaced in <f>. We then project these vector 
velocities onto the line-of-sight at each pixel in a 
4096-by-4096 image of the solar disk as seen from 
1 AU. The average Doppler velocity in each pixel 
is calculated by finding the line-of-sight projection 
at a 7-by-7 array of sub-pixels using bi-cubic in¬ 
terpolation from points on the vector velocity grid 
and then averaging over those sub-pixels that are 
on the visible disk. 

The line-of-sight velocity at a point ( 9 , </>) is 
given by 


Vi os (9, <t>) —V r (9 , <j>) [sin B 0 cos 9 + cos B 0 sin 9 cos (j)\ + 
Vs(9,<j>) [sin B 0 sin 9 — cos B 0 cos 9 cos (j>\ + 
<t>) [cos Bo sin 0] 

( 12 ) 


where Bq is the latitude at disk center (or equiv¬ 
alently the tilt of the Sun’s north pole toward the 
observer) and velocities away from the observer 
are taken to be positive. 

We also add, at each pixel, a convective blue 
shift signal that is a function of the heliocentric 
angle from disk center. This signal complicates the 
extraction of the Doppler signal due to the merid¬ 
ional flow. If these flows are not properly char¬ 
acterized and removed, they leave behind large- 
scale Doppler signals that impact the spectrum 
at low wavenumbers. By also including this “fic- 
tious” signal we faithfully reproduce the full data 
and analysis associated with the HMI data. 



Spatial Wavenumber (pixel ') 


Fig. 6.- The modulation transfer function 
(MTF) for the HMI instrument (red line) and the 
functional fit used for the simulation (black line). 


A point-spread-function (PSF) was developed 
to remove stray light from HMI data and is used 
here in the data simulation. First, a modula¬ 
tion transfer function (MTF) was obtained from 
ground-based calibra tion data using field stops 
( Wachter et alJ[201 2) in order to describe the pre¬ 
launch instrumental optics. We fit the MTF 
with an exponential convolved with the ideal op¬ 
tical transfer function (OTF), often known as a 
‘chat’ function. The exponential and ‘chat’ func¬ 
tion in the frequency domain are equivalent to a 
Lorentzian and Airy function in the spatial do¬ 
main. Post-launch data including solar aureole, 
lunar eclipse and Venus transit events were used 
to evaluate how well the PSF was able to repro¬ 
duce the observed scattering. 

The final value of gamma (the wavenumber ex- 
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ponent) in the exponential function was deter¬ 
mined through least squares fitting of the transit 
of Venus data on 2012.06.05 from the side cam¬ 
era data. The final, non-ideal form of the MTF 
and subsequent PSF was further modified for the 
following reasons. Lunar eclipse images, where a 
large portion of the solar image is obscured by 
the lunar disk, enabled a measure of the large- 
scale scattering. One such eclipse occurred on 
2010.10.07 and the HMI continuum intensity fil- 
tergram shows a light level of 0.34% of the disk- 
center continuum intensity for a position 200 pix¬ 
els onto the lunar disk. We found that the light 
level tended towards a constant far away instead 
of continually decreasing with increasing distance 
from the solar limb. This motivated an additional 
term to be added to the PSF to fit the tail of the 
distribution. If we only considered light scattered 
from 10 ,/ away then the additional term would not 
be necessary, meaning that using only the Venus 
transit data was not sufficient in order to emu¬ 
late the large-scale scattering/dust on the optics. 
(Note that the PSF we developed differs from pre¬ 
vious stray light removal efforts since we do not 
use a single Gaussian or sum of Gaussians as the 
central mathematical component.) 

Deconvolved images were compared to the orig¬ 
inals and it was found that the minimum inten¬ 
sity of a sunspot umbrae decreased from being 
5.5% of the nearby quiet-Sun continuum inten¬ 
sity in the original image to being 3.3% in the 
deconvolved image while the granulation contrast 
doubled, with a standard deviation of the inten¬ 
sity in the quiet-Sun being 3.7% of the average in 
the original image and 7.2% of the average in the 
deconvolved image. 

In our data simulation we average the line-of- 
sight velocity over all sub-pixels within a pixel that 
fall on the visible disk of the Sun and then con¬ 
volve that image with the PSF found for the HMI 
instrument. This convolution is implemented by 
multiplying the FFT of the image with the FFT 
of the PSF (the MTF) and then performing an in¬ 
verse FFT on the result. The MTF found for HMI 
is shown in Figure [6] along with the functional fit 
(which is what we actually use) given by 


MTF(k) = exp(—fc 5 )(l — fc/2) 3 sin(fc7r/2)/(fc7r/2) 

(13) 


An example of a Doppler image from the data 
simulation after the signals due to the axisymmet- 
ric flows and the convective blue shift are removed 
is shown in Figure [7) 



Fig. 7. — Simulated Doppler velocity image. The 
line-of-sight velocity signal ranges from -1000 m 
s^ 1 (dark blue) to +1000 m s _1 (dark red). 


4. DOPPLER VELOCITY SPECTRUM 


We calculate Dopple r velocity spectra following 
Hathaway et al.l ( 2000h . We remove active regions 


in the HMI data by locating pixels where the mag¬ 
netic flux density exceeds 1200 G and then mask¬ 
ing out those pixels along with all adjacent pix¬ 
els within a 10 pixel radius. The full disk HMI 
data, its mask, and the simulated data are pro¬ 
jected onto maps in heliographic longitude and lat¬ 
itude. The mask edges are smoothed by convolv¬ 
ing the mask with a cosine bell having a radius of 
10 pixels. The masked data is then projected onto 
spherical harmonics by taking Fourier transforms 
in longitude and Legendre transforms in latitude 
to obtain complex spectral coefficients, A that 
reproduce the masked data when multiplied by the 
spherical harmonics and summed over t and m. 


The HMI data are comprised of 59 Doppler im¬ 
ages acquired during November 2012 at the hourly 
intervals when the spacecraft motion toward or 
away from the Sun was less than 300 m s -1 . The 
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simulated data are comprised of 18 Doppler images 
constructed with different sets of random numbers 
for the spectral phases. The average velocity spec¬ 
tra are shown in Figure [8] 



Fig. 8.- The Doppler velocity spectra for the 
simulated data (black line) and the HMI data (red 
line). The dotted red line shows the spectrum ob¬ 
tained from the HMI data when the artifacts are 
not removed. 

These spectral values are calculated by taking 
the square root of the wavenumber times the ve¬ 
locity power per wavenumber: 


V(t) 


N 


i 

i y m 2 

m=—£ 


(14) 


The Doppler velocity spectrum rises to a peak 
at £ ~ 120 (consistent with typical supergranules 
with diameters of ~ 35 Mm), drops slightly and 
then rises again to a second peak at t ~ 3500 
(consistent with granules with diameters of ~ 1.2 
Mm). 

Note that this spectrum is for the Doppler ve¬ 
locities. It includes foreshortening at the limb, 
which limits the signal at high wavenumbers, and 
it includes line-of-sight effects, which exclude hor¬ 
izontal velocities at disk center and only include 
one horizontal velocity component across the rest 
of the disk. In spite of these caveats, it is nonethe¬ 
less clear that the spectrum is continuous with 
only two primary features — peaks representative 
of granules and supergranules. 

While cells of all sizes are included, there isn’t 
a spectral feature that indicates a preference for 


mesogranules of any particular size. The bottom 
of the spectral dip at l ~ 400 represents cells with 
diameters of ~ 10 Mm — a quoted size for meso¬ 
granules. 

We have forced the simulation spectrum to drop 
smoothly at low wavenumbers. While this leaves 
an apparent excess of power in the HMI data at 
wavenumbers below 30, we attribute this signal 
excess to a less than perfect removal of the image 
artifacts. This is supported by the low wavenum¬ 
ber feature seen in the data when the artifacts are 
not removed (dotted red line in Figure©- 


5. RADIAL FLOW ANALYSIS 


The radial flow component represented by R™ 
can be separated from the horizont al flow com po¬ 
nents using the method described bv lHathawav et al 
( 2002h . The Doppler velocity due to the horizon¬ 
tal flows drops to zero at disk center, where all 
horizontal velocities are transverse to the line-of- 
sight, while the radial velocity has its strongest 
contribution to the Doppler signal at disk center. 

The line-of-sight velocity is given by 


V los (e,<j>) =V r {9,^)cosp + V h i(9,<j))smp (15) 

where Vhi is the component of the horizontal flow 
velocity aligned with the radius vector from disk 
center, and p is the heliocentric angle from disk 
center. The orthogonal component of the horizon¬ 
tal flow is transverse to the line-of-sight and does 
not contribute to the Doppler signal. Squaring the 
Doppler velocity and averaging it in annuli about 
disk center at a given p gives 


V L (P) = V r cos 2 p + V% x sin 2 p+ 

2V r Vhi cos p sin p 

where the overbars represent the averages. We 
expect the last term in this equation to be neg¬ 
ligible. There is no physical or geometric reason 
that V r Vhi should have a nonzero average on an 
annulus. For example, the up flow at the center of 
a cell should produce horizontal flows away from 
disk center on one side and toward disk center on 
the other - giving a net zero average. Dropping 
this term, we rewrite the equation as 














V 1 2 0S (P) = K 2 + [VZ 1 -V?] s m>p (17) 

This suggests a method for separating the radial 
flow from the horizontal flow component in the 
Doppler velocity signal. If we plot the average 
Doppler velocity squared as a function of sin 2 p, 
then the y-axis intercept gives the radial velocity 
squared and the slope gives the difference between 
the square of the horizontal component and the 
square of the radial component. 

We need to know the contributions to the 
Doppler signal from the radial and horizontal 
flows as functions of wavenumber. To do so we 
reconstruct the Doppler signal from the spherical 
harmonic spectral coefficients, A™, multiplied by 
spectral filters to isolate Doppler features of dif¬ 
ferent sizes. We chose a set of Gaussian filters of 
width A£ = 64 centered on multiples, n, of £ = 128 
from t = 0 to t = 4096 with 

Filter n (t) = exp[-(£ - nl28) 2 /(2A£ 2 )] (18) 

The average Doppler velocity squared for each 
of these 33 filtered Doppler images is fit to a second 
order polynomial in sin 2 p (foreshortening away 
from disk center introduces curvature to these 
data). The y-axis intercept yields the radial flow 
velocity squared while the slope yields the differ¬ 
ence between the squared radial and squared hori¬ 
zontal flow velocities. The radial flow components 
are shown in Figure [9] as functions of t for both 
the HMI and the simulated data. The ratios of 
the radial flow component to the horizontal flow 
component are shown in Figure HU] as functions of 
£ for both the HMI and the simulated data. 

Note that this analysis compares the radial flow 
velocity to just one component of the horizon¬ 
tal flow velocity. While it is clear that the ra¬ 
dial flow is very weak at low wavenumbers and 
becomes comparable to the horizontal flows at 
large wavenumbers, a quantitative value is best 
obtained from the vector velocities in the simula¬ 
tion. 

The simulation has radial flow velocities that 
are just 3% of the total flow velocity at the lowest 
wavenumbers (£ ~ 10). This ratio rises to 6% at 
supergranule wavenumbers (l ~ 120) and to 45% 
at granule wavenumbers (£ ~ 4000). 



Fig. 9.— The radial flow component of the 
Doppler velocity for the simulated data (black 
line) and the HMI data (red line). 



Fig. 10. The ratio of the radial flow component 
to the horizontal flow component of the Doppler 
velocity for the simulated data (black line) and the 
HMI data (red line). 

6. TOROIDAL/POLOIDAL FLOW ANAL¬ 
YSIS 

We have devised a new analysis technique to 
separate the toroidal flow component from the 
poloidal flow component. We construct two ve¬ 
locity gradient images from the prepared data by 
calculating the pixel-wise gradient of the Doppler 
velocity (a scalar) at each point in the image and 
finding the projection of that gradient in the disk- 
radial direction (from disk center toward the limb) 
to produce one image (the Grad p image) and the 
projection of the gradient in the tangential direc¬ 
tion (clockwise about disk center) to produce the 
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second image (the Grad© image). 

A poloidal flow away from a convection cell cen¬ 
ter will contribute to the Grad p image but not to 
the Grad© image along lines that pass radially 
and tangentially through the cell center. Like¬ 
wise, a toroidal flow around the center of a vor¬ 
tex anywhere on the disk will contribute to the 
Grad© image but not to the Grad p image along 
these same lines. Unfortunately, a radial flow will 
contribute to both images and both poloidal and 
toroidal flows will contribute to both images away 
from the lines that pass through cell centers in the 
radial and tangential directions. While the separa¬ 
tion into poloidal and toroidal components is not 
clean and simple, the data simulation allows us 
to nonetheless quantify the contributions of these 
flows to the Grad p and Grad© images. 

We produce these Grad p and Grad© images for 
both the HMI and the simulated data and then 
find their spherical harmonic spectral coefficients 
using the same procedure used for the Doppler 
velocity spectrum in Section 4. 



Fig. 11. The spectra from the Grad p im- 
ages(solid lines) and Grad© images (dotted lines) 
for the simulated data with both poloidal and 
toroidal flows (black lines), the HMI data (red 
lines), and the simulated data in which the hori¬ 
zontal flows are purely poloidal without a toroidal 
component (green lines). 

The results are shown in Figure [TT| We il¬ 
lustrate the sensitivity of the method to poloidal 
and toroidal flows using two simulations. Both 
simulations have horizontal flows that produce 
the spectral match shown in Figure [8] and radial 
flows that produce the spectral matches shown 


in Figures [9] and 1101 One simulation, shown in 
black in Figure [TT] has both poloidal and toroidal 
flows adjusted to match the HMI Grad p and Grad© 
spectra (shown in red). The second simulation 
(shown in green), has only poloidal flows without 
a toroidal flow component. The excess velocity 
at all wavenumbers in the Grad p image spectrum 
from this simulation indicates that a substantial 
fraction of the observed horizontal flows must be 
toroidal. The deficient velocity at wavenumbers 
below i ~ 30 in the Grad© image spectrum from 
this simulation suggests that the toroidal flows 
dominate at these lower wavenumbers. 

7. CONCLUSIONS 

As a result of these analyses, we have simulated 
flows that closely reproduce the analysis results 
from the HMI data seen in the spectral matches 
in the previous sections. The simulations require 
prescribing all three spectral coefficients (radial, 
poloidal, and toroidal) for all wavenumbers from 1 
to 4096. The total velocity spectrum and its three 
components are plotted in Figure [T2] as functions 
of spherical harmonic degree l. Note that these 
spectra are from the vector velocities over the en¬ 
tire solar surface and that they are not impacted 
by the foreshortening or line-of-sight projection ef¬ 
fects that impact the spectra from the Doppler 
images. 



Fig. 12.— The total velocity spectrum (black line) 
from the data simulation. The contribution from 
the poloidal flow is shown in red. The contribution 
from the toroidal flow is shown in blue. The con¬ 
tribution from the radial flow is shown in green. 

As with the Doppler Velocity spectrum, we see 
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that the total velocity spectrum rises to the peak 
at i ~ 120, consistent with supergranules, drops 
slightly at higher wavenumbers and then rises to 
the second peak at £ ~ 3500, consistent with 
granules. The two velocity components coupled 
by the mass continuity equation, the radial and 
poloidal components, have pronounced peaks con¬ 
sistent with supergranules while the toroidal com¬ 
ponent does not. 

None of the components show any spectral fea¬ 
tures indicating a preferred cell size at the lowest 
wavenumbers representing giant cells or at the in¬ 
termediate wavenumbers representing mesogran- 
ules. 


We conclude that, in general, the convection 
spectrum is a rising spectrum with velocities in¬ 
creasing linearly (slope of ~ 1 on this log-log 
plot) with wavenumber t to an ultimate peak at 
wavenumbers of l ~ 3500 consistent with the pho- 
tospheric granulation. However, this linear rise is 
punctuated by a peak at t ~ 120 consistent with 
supergranules. 


The origin of this su pergranule peak is still un¬ 
certain. Simon fc Leighton! (1964) suggested that 
cells the size of supergranules are driven by the 
changes in mean molecular weight at depths where 
helium becomes ionized (depths of ~ 7 Mm for the 
first ionization of He and ~ 30 Mm for th e second 
ionization of He). [November et al.1 ( 1981 1 isolated 
mesogranule-sized cells and suggested that meso- 
granules were driven by the first ionization of He 
and s uperg ranules by the second ionization of He. 
iRastl ( 2003 1 suggested that cells the size of meso- 
granules and supergranules would be driven by the 
collective interactions between downdrafts at cell 
boundaries. However, radiative-MHD simulations 
of co nvectio n in t he outer 20 Mm of the Sun by 
Stein et al. ( 2011 1 include all of these processes, 


but do not yield any dominant cells sizes in the 
photosphere other than granules. (These simula¬ 
tion do, however, indicate that the smaller struc¬ 
tures disappear in deeper layers, leaving velocity 
structures with sizes proportional to the depth.) 

The relative strength of the radial flow velocity 
to the total flow velocity increases from 3% at the 
lowest wavenumbers to nearly 50% at the highest 
wavenumbers. This is in rough agreement with the 
re sults found with t he low er resolution MDI data 
bv lHathawav et ah ( 2002h . but extends the results 
to much higher wavenumbers, and slightly lowers 


the radial flow velocities at the lowest wavenum¬ 
bers (where imaging artifacts were removed from 
the HMI data). 

The toroidal flow component dominates at low 
wavenumbers and becomes comparable to the 
poloidal component at i ~ 30, where velocities 
of ~ 70 m s -1 and cell sizes of ~ 150 Mm give a 
characteristic time scale of ~ 24 days — compa¬ 
rable to the Sun’s rotation period. 

The spectrum shown in Figure [12] is derived 
from Doppler velocity data obtained with HMI. 
This spectrum can be compared to spectra ob¬ 
tained by other means. 



Fig. 13. The total velocity spectrum (red line) 
from the data simulation is plotted for comparison 
with the spect rum from the vector ve locities found 
for giant cells ( Hathaway et al .1201 .'ll) and the spec¬ 
trum from the vector velocities in the StellarBox 
simulation of solar granulation. 


Hathaway et al.l ( 2013il measured the giant cell 


sized flows by a local correlation tracking method 
on the Doppler velocity features (largely super¬ 
granules). This analysis produced maps of the 
horizontal vector velocities for all solar longitudes 
and most solar latitudes (the extreme polar re¬ 
gions were inaccessible). The velocity spectrum 
averaged over 34 solar rotations from May 2010 
through November 2012 is shown in Figure [13] by 
the black line at low wavenumbers. The giant cell 
flow velocities are somewhat smaller that those in 
our simulation, but the slope agrees very well. Our 
data simulation may overestimate the velocity at 
these low wavenumbers by having to fit to HMI 
data that is still impacted by the presence of imag¬ 
ing artifacts. On the other hand, the giant cell flow 
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velocities found bv lHathawav et ah ( 2013 ) may be 
underestimated due to averaging the moving and 
evolving flows at each point over the 10-13 days 
that point is visible in HMI data. 

Numerical simulations for the three-dimensional 
radiative magnetohydrodynamics have progressed 
to the point that they now reproduce the appear¬ 
ance, the time evolution, and the emergent radi ¬ 


ation associated with granules (IKitiashvili et al 


20121 ). We have compared our Doppler veloc¬ 


ity spectrum with the spectrum calculated from 
the three-dimensional radiative hydrodynamics 
simulations obtain ed using the StellarBox code 
I Wrav et al. 20151 ). which takes into account the 


effects of ionization, sub-grid-scale turbulence, in¬ 
ternal structure and chemical composition. For 
analysis we used a fully developed hydrodynamic 
run in a slab 80 by 80 Mm in the horizontal, 25 Mm 
in depth below the photosphere, and 1 Mm of the 
atmospheric layer with the spatial grid resolution 
of 100 km. The photospheric velocity spectrum 
from this simulation (shown in Figure [13] by the 
blue line at high wavenumbers) is in very good 
agreement with the HMI Doppler velocity spec¬ 
trum in the range of the spherical harmonic degree 
from 200 to 4000. However, the supergranulation 
bump is missing in these nu merica l simul ation s (as 
it was in the simulations of Stein et ahl ( 201lh ). 

Understanding of the origin of supergranu¬ 
lation is still a significant challenge in solar 
physics. Origi ns based o n the ionization of helium 
( Simon fc Leightonll964l) or the collective interac¬ 
tions between downdrafts ( Rast 2003) now seem 
unlikely based on numerical simulations. The su¬ 
pergranules are imbedded in the Sun’s surface 
shear layer (which extends to a depth of ~ 50 
Mm) and are characterized by the magnetic net¬ 
work that forms around the periphery of each cell. 
Perhaps their origin lies with the influence of the 
Sun’s rotation or magnetic field on these convec¬ 
tive flows. 
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